Anomalous parity asymmetry of WMAP 7 year power spectrum data at low 
multipoles: is it cosmological or systematics? 



O 

(N 

Oh- 

(N 



o 

u 

Or 
i < 

o 

to 



> 
00 



(N 
O 

o 



X 



J aiseung Kirr0 and Pavel Naselsky 
Niels Bohr Institute & Discovery Center, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 

(Dated: September 3, 2010) 

It is natural to assume a parity neutral Universe, and accordingly no particular parity preference 
in CMB sky. However, our investigation based on the WMAP 7 year power spectrum shows there 
exist large-scale odd-parity preference with high statistical significance. We also find the odd-parity 
preference in WMAP7 data is slightly higher than earlier releases. We have investigated possible 
origins, and ruled out various non-cosmological origins. We also find primordial origin requires 
Re[$(k)]| <C |Im[$(k)]| for k < 22/rjo, where rjo is the present conformal time. In other words, it 
requires translational invariance in primordial Universe to be violated on the scales larger than 4 Gpc. 
The Planck surveyor, which possesses wide frequency coverage and systematics distinct from the 
WMAP, may allow us to resolve the mystery of the anomalous odd-parity preference. Furthermore, 
polarization maps of large sky coverage will reduce degeneracy in cosmological origins. 

PACS numbers: 95.85.Sz, 98.70.Vc, 98.80.Cq, 98.80.Es, 98.80.-k 



I. INTRODUCTION 

For the past years, there have been great successes in 
measurement of CMB anisotropy by ground and satel- 
lite observations [3-Hlj. CMB anisotropy, which are as- 
sociated with the inhomogeneity of the last scattering 
surface, provides the deepest survey so far and allows 
us to constrain cosmological models significantly. Since 
the initial release of WMAP data, the WMAP CMB 
sky data have undergone scrutiny, and various anoma- 
lies have been found and reported [12-33]. In particular, 
CMB anisotropy at low multipoles are associated with 
scales far beyond any existing astrophysical survey, and 
therefore CMB anomalies at low multipoles may hint new 
physical laws at unexplored large scales. 

CMB sky map may be considered as the sum of even 
and odd parity functions. Previously, Land and et al. 
have noted the odd point-parity preference of WMAP 
CMB data, but found its statistical significance is not 
high enough [24| . In our previous work, we have ap- 
plied a slight different estimator to WMAP 5 year power 
spectrum data, and found significant odd point-parity 
preference at low multipoles [31] . In this paper, we in- 
vestigate the recently released WMAP 7 year data up 
to higher multipoles, and discuss origins of the observed 
odd-parity preference. 

This paper is organized as follows. In Section [Til 
we discuss the anomalous odd-parity preference of the 
WMAP data. In Section IIIIl we implement cosmologi- 
cal fitting by excluding even or odd low multipole data. 
In Section IIVI and [V] we investigate non-cosmological 
and cosmological origin. In Section IVI1 we summarize 
our investigation and discuss prospect. In Appendix (X[ 
we briefly review statistical properties of Gaussian CMB 
anisotropy. 



II. PARITY ASYMMETRY OF THE WMAP 
DATA 

CMB anisotropy sky map may be considered as the 
sum of even and odd parity functions: 



where 



T(n) = T+(n)+r-(n), 
T(n)+T(-n) 



T+(n) 

r-(n) 



T(n)-T(-n) 



(1) 

(2) 
(3) 



Taking into account the parity property of spherical har- 
monics Y/ m (n) = (—1)' Yi m (— n) [3J], we may easily show 



T+(n) = J2 ai m Y lm (n), 

l—2n,m 

T~(n) = 2J aim,Yi m (n), 



(4) 
(5) 



i=2n-l.' 
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where n is an integer. Therefore, significant power asym- 
metry between even and odd multipoles indicates a pref- 
erence for a particular parity. Hereafter, we will denote 
a preference for particular parity by 'parity asymmetry'. 
In Fig. [1] we show the WMAP 7 year, 5 year, 3 yea r 
data and the WMAP concordance model 0, H 
We may see from Fig. [1] that the power spectrum of 
WMAP data at even multipoles tend to be lower than 
those at neighboring odd multipoles. In Fig. [21 we show 
+ 1)/2tt (Q wmap - C ACDM ) for low multipoles. 
Since we expect random scattering of data points around 
a theoretical model, we expect the distribution of dots in 
Fig. [2] to be symmetric around zero. However, there are 
only 5 points of positive values among 22 points in the 
case of WMAP7 or WMAP 5 data. Therefore, we may 
see that there is the tendency of power deficit (excess) at 
even (odd) multipoles, compared with the ACDM model. 
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FIG. 1: CMB power spectrum: WMAP 7 year data (blue), 
WMAP 5 year data (green) and WMAP 3 year data (red), 
ACDM model (cyan) 
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2: ( — l)'x difference between WMAP power spectrum 
and ACDM model 



Taking into account 1(1 + 1)C; 
the following quantities: 



const, we may consider 



1=2 

(max 

P- = £ 2^(1 -(-I)') J(l + l)/27rC,, (7) 

1=2 

where P + and P~ are the sum of 1(1 + l)/2w Ci for even 
and odd multipoles respectively. Therefore, the ratio 
P + /P~ is associated with the degree of the parity asym- 
metry, where the lower value of P + / P~ indicates odd- 
parity preference, and vice versa. It is worth to note that 
our estimator of the parity asymmetry does not possess 
explicit dependence on the underlying theoretical model 



in the sense that the term (7 ; WMAP 
our estimator. 



C, ACDM is absent in 



FIG. 3: P+ /P~ of WMAP data and ACDM 



In Fig. [3j we show the P+ /P~ of WMAP data, and 
a ACDM model for various / max - As shown in Fig. G2 
P + /P~ of WMAP data are far below theoretical val- 
ues. Though the discrepancy is largest at lowest Z max , 
its statistical significance is not necessarily high for low 
I, due to large statistical fluctuation. In order to make 
rigorous assessment on its statistical significance at low 
I, we are going to compare P + /P~ of WMAP data 
with simulation. We have produced 10 4 simulated CMB 
maps HEALPix Nside=8 and Nside=512 respectively, via 
map synthesis with a; TO randomly drawn from Gaussian 
ACDM model. We have degraded the WMAP processing 
mask (Nside=16) to Nside=8, and set pixels to zero, if 
any of their daughter pixels is zero. After applying the 
mask, we have estimated power spectrum 2 < I < 23 
from simulated cut-sky maps (Nside=8) by a pixel-based 
maximum likelihood method P, HI [H. At the same 
time, we have applied the WMAP team's KQ85 mask to 
the simulated maps (Nside=512), and estimated power 
spectrum 2 < I < 1024 by pseudo C; method (iol. Hljj. 
In the simulation, we have neglected instrument noise, 
since the signal-to-noise ratio of the WMAP data is quite 
high at multipoles of interest (i.e. I < 100) [l|, [2|. Us- 
ing the low I estimation by pixel-maximum likelihood 
method and high I estimation by pseudo C/ method, we 
have computed P + / P~ respectively for various multi- 
pole ranges 2 < I < l max , and compared P + / P~ of the 
WMAP data with simulation. In Fig. |4l we show p- value 
of WMAP7, WMAP 5 and WMAP3 respectively for var- 
ious / m ax, where p- value denotes fractions of simulations 
as low as P + /P~ of the WMAP data. As shown in Fig. 
IH the parity asymmetry of WMAP7 data at multipoles 
(2 < I < 22) is most anomalous, where p- value is 0.0031. 
As shown in Fig. SJ the statistical significance of the par- 
ity asymmetry (i.e. low p-value) is getting higher, when 
we include higher multipoles up to 22. Therefore, we may 
not attribute the odd parity preference simply to the low 
quadrupole power, and find it rather likely that the low 
quadrupole power is not an isolated anomaly, but shares 
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FIG. 4: Probability of getting P + /P~ as low as WMAP data 
for multipole range 2 < I < Z max - 



an origin with the odd parity preference. 

TABLE I: the parity asymmetry of WMAP data (2 < I < 22) 



data 


P +/p- 


p- value 


WMAP7 


0.7076 


0.0031 


WMAP5 


0.7174 


0.0039 


WMAP3 


0.7426 


0.0061 



In Table HI we summarize P + / P and p- values of 
WMAP7, WMAP5 and WMAP 3 for l max = 22. As 
shown in Fig. [4] and Table U the odd-parity preference 
of WMAP7 is most anomalous, while WMAP7 data are 
believed to have more accurate calibration and less fore- 
ground contamination than earlier releases. QJ-0, l42j . 
In Fig. we show cumulative distribution of P + / P~ 
for 10 simulated maps. The values corresponding to 
P + /P~ of WMAP data are marked as dots. We have 
also compared P + / P~ of the WMAP7 with whole-sky 
simulation (i.e. no mask), and obtained p-value 0.002 
for Z max = 22. The lower p- value from whole-sky simula- 
tion is attributed to the fact that statistical fluctuation 
in whole-sky Ci estimation is smaller than that of cut-sky 
estimation. 

In the absence of strong theoretical grounds for the 
parity asymmetry (2 < I < 22), we have to take into 
account our posteriori choice on Z max , which might have 
enhanced the statistical significance. In order to do that, 
we have produced whole-sky Monte-Carlo simulations, 
and retained only simulations, whose p-value is lowest at 
lmax = 22. The p- values have been estimated by com- 
paring them with simulations, which are produced sepa- 
rately. Once we have retained 10 4 simulations, we have 
compared them with WMAP7 data, and found that only 
fraction 0.0197 of retained simulations have P + / P~ as 
low as WMAP7 data. The statistical significance of the 
parity asymmetry (2 < I < 22) is reduced substantially 
by accounting for the posteriori choice on l max . However, 
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FIG. 5: Parity asymmetry at multipoles (2 < I < 22): cumu- 
lative distribution of P + / P~ for 10 4 simulated maps (cyan), 
P+/P~ of WMAP 7 (blue), WMAP 5 (green) and WMAP 3 
(red) 

it is still significant. 

We have also investigated the parity asymmetry with 
respect to mirror reflection (i.e. I + m=cven or odd) re- 
spectively in Galactic coordinate and Ecliptic coordinate. 
However, we find the statistical significance is not as high 
as the point-parity asymmetry (i.e. Z=even or odd). 



III. THE POWER CONTRAST AND THE 
ACDM MODEL FITTING 

The parity asymmetry discussed in the previous sec- 
tion is explicitly associated with the angular power spec- 
trum data, which are used extensively to fit cosmolog- 
ical models. Noting the significant power contrast be- 
tween even and odd low multipoles, we have investigated 
cosmological models respectively by excluding even or 
odd low multipole data (2 < I < 22) from total dataset. 
Hereafter, we denote CMB data of even(odd) multipoles 
(2 < I < 22) plus all high multipoles as 'ZV('Ai') re- 
spectively. For a cosmological model, we have consid- 
ered ACDM + SZ effect + weak-lensing, where cosmo- 
logical parameters are A £ {fib, fl c , t, n s , A s , A sz , Hq}. 
For data constraints, we have used the WMAP 7 year 
power spectrum data [![. In order to exclude even or 
odd low multipoles, we have made slight modifications 
to the likelihood code provided by the WMAP team, 
and run a CosmoMC with the modification, [l], 0, [lij . 
In Fig. [6J we show the marginalized likelihoods of pa- 
rameters. As shown in Fig. |5l the parameter likelihood 
imposed by D3 seem to differ from those of others. In 
Table |H] we show the bestfit parameters and 1 a con- 
fidence intervals, where A, A2 and A3 denote the best- 
fit values of the full data, D2 and D% respectively. As 
mentioned above, there exist some level of tension be- 
tween D3 and the full data. However, the Bayes Factor 
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FIG. 6: Marginalized likelihood of cosmological parameters: 
results are obtained respectively with/without even or odd 
multipole data (2 < I < 22). 



TABLE II: cosmological parameters (ACDM ± sz ± lens) 



A 



A 2 



A 



3 



lO! 



n c h 2 

T 

n s 
l[W 10 A s 
H 

As?, 



0.0226 ± 0.0006 0.0228 ± 0.0006 0.0221 ± 0.0006 



0.112 ±0.006 
0.0837 ±0.0147 
0.964 ± 0.014 
3.185 ±0.047 
70.53 ± 2.48 
1.891^ 



0.11 ±0.006 
0.0879 ± 0.015 
0.974 ±0.015 
3.165 ±0.049 

71.43 ±2.51 



I +0.109 
-1.891 



1.469 



+0.541 



0.116 ±0.006 
0.087 ± 0.0147 

0.95 ±0.015 
3.246 ± 0.048 

68.07 ±2.53 

1.558:! 



2+0.442 
-1.558 



C(X\D 3 )/C(X 3 \D 3 ) = exp(-3733.124± 3732.791) = 0.71 
shows that it is 'not worth than a bare mention', accord- 
ing to Jeffreys' scale (45[. 



IV. NON-COSMOLOGICAL ORIGINS 

In this section, we are going to investigate non- 
cosmological origins such as asymmetric beams, instru- 
ment noise, foreground and cut-sky effect. 



FIG. 7: the parity asymmetry in the presence of beam asym- 
metry: The dots denotes (P + ,P~) of CMB maps simulated 
with asymmetric beams. The dashed lines are plotted with 
slopes corresponding to the P + /P~ of ACDM model (red) 
and WMAP7 data (green) respectively. The alphanumeric 
values at the lower right corner denote the frequency band 
and D/A channel. 



sumed symmetric beams in the power spectrum estima- 
tion [HUES, EH- We have investigated the association of 
beam asymmetry with the anomaly, by relying on simu- 
lated maps provided by (47j . The authors have produced 
10 simulated maps for each frequency and Differencing 
Assembly (D/A) channels, where the detailed shape of 
the WMAP beams and the WMAP scanning strategy are 
taken into account [471 ]. From simulated maps, we have 
estimated P + and P~, where we have compensated for 
beam smoothing purposely by the WMAP team's beam 
transfer function (i.e. symmetric beams). In Fig. [71 
we show P + and P~ values of the simulated maps, and 
the dashed lines of a slope corresponding to P + / P~ of 
ACDM and WMAP7 data respectively. As shown in Fig. 
[7J we do not observe the odd-parity preference of WMAP 
data in simulated maps. Therefore, we find it hard to at- 
tribute the odd-parity preference to asymmetric beams. 



A. asymmetric beam 

The shape of the WMAP beams are slightly asym- 
metric 0, HI, H3|, while the WMAP team have as- 



Thcrc exist instrument noise in the WMAP data. Es- 
pecially, 1/f noise, when coupled with WMAP scanning 
pattern, may result in less accurate measurement at cer- 
tain low multipoles [H, El: EH . In order to investigate 
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the association of noise with the anomaly, we have pro- 
duced noise maps of WMAP7 data by subtracting one 
Differencing Assembly (D/A) map from another D/A 
data of the same frequency channel. In Fig. [HI we show 
P + and P~ values of the noise maps. As shown in Fig. 
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FIG. 8: the parity asymmetry of the WMAP noise: The dots 
denotes (P + ,P~) of noise maps, and alphanumeric values in 
the legend denote the frequency band and the pair of D/A 
channels used. Two dashed lines are plotted with the slope 
corresponding to P + jP~ of white noise and WMAP7 data 
respectively. 

[HI the noise maps do not show odd-parity preference, 
but their P + /P~ ratios are consistent with that of white 
noise (i.e. Ci — const.). Besides that, the signal-to-noise 
ratio of WMAP temperature data is quite high at low 
multipoles (e.g S/N~ 100 for I = 30) E2, Eg . There- 
fore, we find that instrument noise, including 1/f noise, 
is unlikely to be the cause of the odd-parity preference. 



C. foreground 

There are contamination from galactic and extragalac- 
tic foregrounds. In order to reduce foreground contam- 
ination, the WMAP team have subtracted diffuse fore- 
grounds by template-fitting, and masked the regions that 
cannot be cleaned reliably. For foreground templates 
(dust, free- free emission and synchrotron), the WMAP 
team used dust emission "Model 8", Ha map, and the 
difference between K and Ka band maps 0, l50l - [53T | . In 
Fig. [HI we show the power spectrum of templates. As 
shown in Fig. [HI templates show strong even parity pref- 
erence, which is opposite to that of the WMAP power 
spectrum data. Therefore, one might attribute the odd- 
parity preference of WMAP data to over-subtraction by 
templates. However, we find the error associated with 
template-fitting is unlikely to be the cause of the odd- 
parity preference. Consider spherical harmonic coeffi- 
cients of a foreground-reduced map: 



obs 
a lm 



cmb 
a lm 



ba 



i pi 

Irn ' 




FIG. 9: the power spectra of the templates (synchrotron, Ha, 
dust): plotted with arbitrary normalization. 



where a°^ s , a[?L and baj^ correspond to a foreground- 
cleaned map, a foreground and a template with a fit- 
ting coefficient b. For simplicity, we consider only a sin- 
gle foreground component, but the conclusion is equally 
valid for multi-component foregrounds. Since there is no 
correlation between foregrounds and CMB, the observed 
power spectrum is given by: 



G 



obs 



G 



cmb 



Hm 



I. tpl 

ba hn 



(9) 



As shown Eq. [HI the parity preference should follow that 
of templates because of the second term, provided tem- 
plates are good tracers of foregrounds (i.e. 



const). Nevertheless, Eq. [9] may make a bad approx- 
imation for lowest multipoles, because the cross term 
£ m R« b (a[f n - baft)*] may not be negligible. Be- 
sides that, our argument and the template- fitting method 
itself fail, if templates are not good tracers of foregrounds. 
In order to investigate these issues, we have resorted to 
simulation in combination with WMAP data. Noting the 
WMAP power spectrum is estimated from foreground- 
reduced V and W band maps, we have produced simu- 
lated maps as follows: 



T(n) = T cmb (n) + (V(h) - W(h))/2, 



(10) 



(8) 



where V(h) and W(h) are foreground-reduced V and W 
band maps of WMAP data. Note that the second term 
on the right hand side contains only residual foregrounds 
at V and W band, because the difference of distinct fre- 
quency channels are mainly residual foregrounds at low I. 
Just as cut-sky simulation described in Sec. [Hi we have 
applied a foreground mask to the simulated maps, and 
estimated the power spectrum from cut-sky by a pixel- 
based maximum likelihood method. In Fig. 1101 we show 
P + and P~ values estimated from simulations. For com- 
parison, we have included simulations without residual 
foregrounds, and dashed lines of a slope corresponding to 
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FIG. 10: the parity asymmetry in the presence of residual 
foregrounds (V— W): Dashed lines are plotted with slopes cor- 
responding to P+/P- of P+/P~ of ACDM (cyan), WMAP7 
data (red). 

P+ /P~ of ACDM model and WMAP7 data. As shown 
in Fig. [TOl the P + / P~ of simulations in the presence of 
residual foregrounds do not show anomalous odd-parity 
preference of WMAP data. Considering Eq. [9] and sim- 
ulations, we find it difficult to attribute the odd-parity 
preference to residual foreground. 

There also exist contamination from unresolved extra- 
galactic point sources [35j]. However, point sources follow 
Poisson distribution with little departure [54| , and there- 
fore are unlikely to possess odd-parity preference. Be- 
sides that, point sources at WMAP frequencies are sub- 
dominant on large angular scales (low I) (35l. IHol IHil. |55| . 

Though we have not find association of foregrounds 
with the anomaly, we do not rule out residual foreground 
completely, because of our limited knowledge on residual 
foregrounds. 
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FIG. 11: Probability of getting P + /P~ as low as the ILC 7 
year, 5 year, and 3 year map at multipole range 2 < I < i max 



P + jP and p- values for Z max = 22. 



TABLE III: the parity asymmetry of WMAP ILC maps (2 < 
I < 22) 



data 


P + /P- 


p- value 


ILC7 


0.7726 


0.0088 


ILC5 


0.7673 


0.0074 


ILC3 


0.7662 


0.0072 



As shown in Fig. [TTJ and Table IIII| we find anoma- 
lous odd-parity preference exits in whole-sky CMB maps 
as well. Therefore, we find it difficult to attribute the 
anomaly to cut-sky effect. 



E. Other known sources of errors 



D. cut sky 

The WMAP team have masked the region that cannot 
be reliably cleaned by template fitting, and estimated 
CMB power spectrum from sky data outside the mask 
[H> 0) US H3- Therefore, we have estimated the p- value 
presented in Sec. [Ill by comparing WMAP data with cut- 
sky simulations. Nonetheless, we have investigated the 
WMAP team's Internal Linear Combination map (ILC) 
map in order to see if the odd-parity preference also ex- 
ists in a whole-sky CMB map. Note that the WMAP 
ILC map provides a reliable estimate of CM B sig nal over 
whole-sky on angular scales larger than 10° [35l 150. l55|. 
We have compared P + /P~ of the ILC maps with whole- 
sky simulations. In Fig. [Til we show p-values of the ILC 
maps respectively for various Z max . As shown in Fig. ITT1 
the odd-parity preference of ILC maps is most anoma- 
lous for Z max = 22 as well. In Table IIII1 we summarize 



Besides contamination discussed in previous sections, 
there are other sources of contamination such as sidelobe 
pickup and so on. In order to investigate these effects, 
we have resorted to simulation produced by the WMAP 
team. According to the WMAP team, time-ordered data 
(TOD) have been simulated with realistic noise, ther- 
mal drifts in instrument gains and baselines, smearing of 
the sky signal due to finite integration time, transmis- 
sion imbalance, and far-sidelobe beam pickup. Using the 
same data pipeline used for real data, the WMAP team 
have processed simulated TOD, and produced maps for 
each differencing assembly and each single year observa- 
tion year. In Fig. [T^l we show the P + and P~ of the 
simulated maps, where the power spectrum estimation is 
made from cut-sky by a pixel-based likelihood method. 
As shown in Fig. [T2l all points are well above P + / P~ of 
WMAP7, and agree with ACDM model. Therefore, we 
do not find definite association of the parity asymmetry 
with known systematics effects. 
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FIG. 12: P+ and P~ of the WMAP team's simulation for V 
and W band data 



V. COSMOLOGICAL ORIGIN 

In this section, we are going to take the WMAP power 
spectrum at face values, and consider cosmological ori- 
gins. Topological models including multi-connected Uni- 
verse and Bianchi VII model have been proposed to 
explain the cold spot or low quadrupole power [56l - l58l |. 
However, the topological models do not produce the par- 
ity asymmetry, though some of them, indeed, predict low 
quadrupole power. Trans-Planckian effects and some in- 
flation models predict oscillatory features in primordial 
power spectrum [59j - [70l j. However, oscillatory or sharp 
features in primordial power spectrum are smeared out 
in translation to the CMB power spectrum ]TlJ - Be- 
sides, reconstruction of primordial power spectrum and 
investigation on features show that primordial power 
spectrum is close to a featureless power-law spectrum 
SISEIllSIziSI. Therefore, we find it difficult to at- 
tribute the anomaly to trans-Planckian effect or extended 
inflation models. We will consider what the odd-parity 
preference imply on primordial perturbation $(k), if pri- 
mordial power spectrum is, indeed, featureless. Using 
Eq. IA21 we may show the decomposition coefficients of 
CMB anisotropy are given by: 



(-»)' 
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(-*) 





OO 7T 
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where we used the reality condition $(— k) = <f>*(k) and 
Yi m (—: n) = (—l) l Yi m (n). Using Eq. [Til it is trivial to 
show, for the odd number multipoles I = 2n — 1, 



air, 



(11) 



OO 7T 



( ^ 1 / dk I ,W k sin« k / ^. k , // (/ 1 -)r,;„(k)Tm[<I>(k) . 





and, for even number multipoles / = 2n, 



air. 



(12) 



OO 7T 



dk / d0 k sin0 k / d</>kffi(fc)r^(k)Re[$(k)]. 



(-1)' 



It should be noted that the above equations are simple 
reformulation of Eq. IA21 and exactly equal to them. 

From Eq. [TTJ and Q21 we may see that the odd-parity 
preference might be produced, provided 



|Re[$(k)]| « |Im[$(k)]| (A; < 22/7,0), 



(13) 



where 770 is the present conformal time. Taking into ac- 
count the reality condition $(— k) = < f > *(k), we may show 
primordial perturbation in real space is given by: 



$(x) 



2 dk d6» k sin6» k 



(14) 




x (Re[$(k)]cos(k-x) 



Im[$(k)] sin(k • x)) . 



Noting Eq. [JJJ] and 1141 we find our primordial Universe 
may possess odd-parity preference on large scales (2/r)o < 
k < 22/770). The odd-parity preference of our primordial 
Universe violates large-scale translational invariance in 
all directions. However, it is not in direct conflict with 
the current data on observable Universe (i.e. WMAP 
CMB data) , though it may seem intriguing. Considering 
Eq. [13] and [JJ] we find this effect will be manifested 
on the scales larger than 2tt rjo/22 « 4Gpc. However, 
it will be difficult to observe such large-scale effects in 
non-CMB observations. If the odd-parity preference is 
indeed cosmological, it indicates we are at a special place 
in the Universe, which may sound bizarre. However, it 
should be noted that the invalidity of the Copernican 
Principle such as our living near the center of void had 
been previously proposed in different context [75|, [7(| . 

Depending on the type of cosmological origins (e.g. 
topology, features in primordial power spectrum and 
Eq. I13[) . distinct anomalies are predicted in polarization 
power spectrum. Therefore, polarization maps of large- 
sky coverage (i.e. low multipoles) will allow us to remove 
degeneracy and figure a cosmological origin, if the parity 
asymmetry is indeed cosmological. 

VI. DISCUSSION 

We have investigated the parity asymmetry of our early 
Universe, using the newly released WMAP 7 year power 
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spectrum data. Our investigation shows anomalous odd- 
parity preference of the WMAP7 data (2 < I < 22) at 3- 
in-1000 level. When we account for our posteriori choice 
on ?max) the statistical significance decreases to 2-in-100 
level, but remains significant. 

There exist several known CMB anomalies at low mul- 
tipole, including the parity asymmetry discussed in this 
paper (1614231 l3ll. 1771. |78| (see [32| for review). We find it 
likely there exist a common underlying origin, whether 
cosmological or not. 

We have investigated non-cosmological origins, and 
ruled out various non-cosmological origins such as asym- 
metric beams, noise and cut-sky effect. We have also 
investigated the WMAP team's simulation, which in- 
cludes all known systematic effects, and do not find defi- 
nite association with known systematics. Among cosmo- 
logical origins, topological models or primordial power 
spectrum of feature might provide theoretical explana- 
tion, though currently available models do not. We also 
find primordial origin requires Re[$(k)]| <C |Im[$(k)] 
for k < 22/?7o, if we consider a simple phenomenologi- 
cally fitting model. In other words, it requires violation 
of translation invariance in primordial Universe on the 
scales larger than 4 Gpc. 

Depending on the type of cosmological origins, distinct 
anomalies are predicted in polarization power spectrum. 
Therefore, we will be able to remove degeneracy in cos- 
mological origins, when polarization data of large sky 
coverage are available. However, at this moment, it is not 
even clear, whether the anomaly is due to unaccounted 
contamination or indeed cosmological. Nonetheless, we 
may be able to resolve the mystery of the large-scale odd- 
parity preference, when data from the Planck surveyor 
are available. 
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Appendix A: Statistical Properties of CMB 
anisotropy 

The CMB temperature anisotropy over a whole-sky is 
conveniently decomposed in terms of spherical harmonics 
Yi m (6, <fi) as follows: 

T(h) = J2aimYlm(n), (Al) 

im 

where ai m is a decomposition coefficient, and n is a sky 
direction. Decomposition coefficients are related to pri- 
mordial perturbation as follows: 

aim = M-*) l J ^m 9 i(k)Y lm (k), (A2) 

where $(k) is primordial perturbation in Fourier space, 
and gi{k) is a radiation transfer function. For a Gaus- 
sian model for primordial perturbation, decomposition 
coefficients satisfy the following statistical properties: 

{ai m ) = 0, (A3) 

{ a lm a l'm') = Cl6wS mm >, (A4) 

where (...) denotes the average over the ensemble of uni- 
verses. Given a standard cosmological model, Sach-Wolf 
plateau is expected at low multipoles [8l|: 1(1 + l)Ci ~ 
const. 
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